Bloch oscillations in Fermi gases 



co 
o 
o 

(N 

wo: 
< 

in 



M. Rodriguez 1 and P. Torma 2 

laboratory of Computational Engineering, P.O.Box 9203, 
FIN-02015 Helsinki University of Technology, Finland 
2 Department of Physics, University of Jyvaskyla, P.O.Box 35, FIN-40014 Jyvaskyla, Finland 

(Dated: February 2, 2008) 

The possibility of Bloch oscillations for a degenerate and superfluid Fermi gas of atoms in an 
optical lattice is considered. For a one-component degenerate gas the oscillations are suppressed 
for high temperatures and band fillings. For a two-component gas, Landau criterion is used for 
specifying the regime where Bloch oscillations of the superfluid may be observed. We show how the 
amplitude of Bloch oscillations varies along the BCS-BEC crossover. 
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The experimental realization of optical lattices for 
bosonic atoms has led to several landmark experiments 
[H S Very recently similar potentials have be- 
come available for trapping the fermionic isotopes as well 
0, . An increase in the superfluid transition temper- 
ature when using potentials created by standing light 
waves has been predicted 0- For trapped cold atoms, 
the famous BCS-BEC crossover problem 0, could be 
studied by tuning the interaction strength between the 
atoms using Feshbach resonances 0,0. In optical lattices 
the whole BCS-BEC crossover could be scanned experi- 
mentally also in an even simpler way by modulating the 
light intensity. We consider Bloch oscillations in these 
systems and show that they can be used as a tool for 
studying the crossover. 

Bloch oscillations are a pure quantum phenomenon oc- 
curing in a periodic potential. They have never been ob- 
served in a natural lattice for electrons as predicted in 
[To| because the scattering time of the electrons by lat- 
tice defects or impurities is much shorter than the Bloch 
period. However, Bloch oscillations have recently been 
observed in semiconductor superlattices 0, for quasi- 
particles penetrating the cores of a vortex lattice in a 
cuprate superconductor 0] , and for periodic optical sys- 
tems such as waveguide arrays [Ijjj. Also cold bosonic 
atoms and superfiuids in optical lattices have been shown 
to be clean and controllable systems well suited for the 
observation of Bloch oscillations 0, 0] • 

Several novel aspects of the physics of Bloch oscilla- 
tions arise for fermionic atoms in optical lattices, i) Im- 
purity scattering can be made negligible, and the parti- 
cle number controlled at will to produce any band filling. 
Even when Bloch oscillations were originally proposed 
for fermions, the effect of the Fermi sea has not played 
a major role. Due to impurity and defect scattering, the 
studies of transport in presence of a constant force have 
focused on drift velocities rather than oscillations. In 
this letter we generalize the semiclassical single particle 
description of Bloch oscillations to arbitrary band fillings, 
ii) The possibility of an oscillating fermionic superfluid 
becomes relevant. We derive the Landau criterion for 



the optical lattice imposing the Cooper pair size to be of 
the order or smaller than the lattice spacing. For solid 
state systems, the Cooper pair radius is usually much 
larger than the lattice spacing and periodicity irrelevant 
for the superfluid, therefore the system is treated as ho- 
mogeneous when calculating supercurrents. We calculate 
the superfluid velocity in the periodic potential and show 
that pairing, leading to smoothening of the Fermi edge, 
suppresses Bloch oscillations. 

Using six counter-propagating laser beams of wave- 
length A, an isotropic 3D simple cubic lattice potential 
can be created which is of the form 
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where Vo is proportional to the laser intensity and a = 
A/2. With the Bloch ansatz the Schrodinger equa- 
tion leads to a band structure in the energy spectrum 
e n (k) . One-component degenerate Fermi gas at low tem- 
peratures can be considered as non-interacting since p- 
wave scattering is negligible and s-wave scattering sup- 
pressed by Fermi statistics. We are interested in high 
enough values of Vq such that tunneling is small and 
tight binding approximation can be applied. The dis- 
persion relation for the lowest band becomes e(k) = 
J[3 — cos(k x a) — cos(k y a) — cos(k z a)], where the band 

width J = -^E R (j^) 3/4 exp (-2^J^\ is obtained 
using the WKB-approximation and Er is the recoil en- 
ergy of the lattice |fj . 

In a two-component Fermi gas, atoms in two different 
hyperfme states (" |, T ") m &y interact with each other. 
The interaction can be assumed pointlike, characterized 
by a scattering length as- The system Hamiltonian 
H = E a Jd 3 Ti>m(T + V)i> a (r) - \g\Jd 3 r^\^ h 
where g — Ai:h 2 as /m can then be mapped to the 



attractive Hubbard model H 



- JE 
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where U = Hir- 

Notc that the BCS (J » U) to BEC (U » J) cross- 
over can be controlled by Vo alone. One-band description 
is used in the Hubbard model also in the case of strong 



interactions [FJ]. We define the limits of the one-band 
approximation for the physical potential Eq.Q by de- 
manding the lowest band gap to be bigger than the effec- 
tive interaction U (note that U > \g\ for the parameters 
of interest). The band gap can be estimated by approxi- 
mating the cosine potential well by a quadratic one. De- 
manding the corresponding harmonic oscillator energy to 

be greater than U gives the condition < (rf~|) ' 
Since a > |ag| imposed by considering on-site interac- 
tions only, the condition is easily valid in general, and for 
the parameters of Fig. [21 in particular. Estimates made 
using exact numerical band gaps in ID support this ar- 
gument. One-band approximation is sufficient because 
larger Vo means steeper optical potential wells which not 
only increase the effective interaction U but also the band 
gaps. 

Bloch oscillations for a single atom can be character- 
ized considering the mean velocity of a particle in a Bloch 
state v(n, k) = (n, k|r|ra, k) given by 

v(n,k) = iv k£n (k). (2) 

When a particle in the Bloch state \n, ko) is adiabati- 
cally affected by a constant external force F = F x x weak 
enough not to induce interband transitions, it evolves 
up to a phase factor into the state |n,k(i)) according 
to k(t) = ko + ^jf 2 . The time evolution has a period 
t~b = h/(\F x \a), corresponding to the time required for 
the quasimomentum to scan the whole Brillouin zone. If 
the force is applied adiabatically, it provides momentum 
to the system but not energy because the effective mass 
(given by m(e) _1 = ^rf^f) is not always positive. For 
optical lattices the force (or tilt: V = F ■ r term in the 
Hamiltonian) can be realized by accelerating the lattice 
0, 0. Using the tight-binding dispersion relation the 
velocity of an atom oscillates like 

v x (t) = sm(k 0x a + F x ta/K). (3) 

For cold bosonic atoms and condensates 0, Q nearly 
all of the population is in the lowest mode of the op- 
tical potential, Eq.(|3J therefore describes the oscillation 
of the whole gas. We generalize the result for the case 
when many momentum states of the band (at T=0, the 
states with wave vector k < kp) are occupied. We calcu- 
late the velocity of the whole gas as the average over the 
normalized temperature-dependent distribution function 
(the Fermi distribution /) of the particles: 

(«*(*)) = /(ko)V fc0x e(k + Ft/ft). (4) 

ko 

Using the tight-binding dispersion relation for the Bloch 
energies we obtain the oscillations shown in Fig. ^ At 
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FIG. 1: The average velocity in 4? units as a function of time 
for a half filled band. The plotted lines correspond to atoms 
in the normal state at different temperatures T — J, 0.5J and 
< 0.1J. Bigger amplitudes correspond to lower temperatures. 

T=0, Eq.Q reduces to 

, , .. Jasm(k xF a) . ( F x ta\ 

Mt)) = T^^- sm {—)- (5) 

This shows that macroscopic coherent oscillation effect 
can still be observed if the band is not full, but the am- 
plitude is suppressed by the band filling k x pa. The effect 
of the temperature can be seen in Fig. ^ the amplitude 
starts to decrease at temperatures of the order T > 0.1J 
but is still non- negligible at half J. These results are valid 
for the one-component degenerate Fermi gas at low tem- 
peratures. In the two-component Fermi gas, atoms in 
the different hyperfine states interact with each other 
which may lead to a superfluid state. Above T c , weak 
interactions can be described by a mean field shift in the 
chemical potential, leading to no qualitative changes in 
Bloch oscillations. Inelastic scattering and consequent 
damping of Bloch oscillations can be described e.g. by 
balance equations |15| . In the following we consider the 
superfluid case where qualitative changes are expected. 

In order to observe Bloch oscillations of the superfluid 
Fermi gas, the critical velocity of the superfluid should 
not be reached before the edge of the Brillouin zone. A 
BCS-superconductor can carry a persistent current q un- 
til a critical velocity, v c = For higher current val- 
ues, even at T — 0, it might be energetically favorable 
to break Cooper pairs and create a pair of quasiparti- 
cles [l^. This costs 2 A in binding energy and decreases 
the Bloch energy by \£ k+q - £ k -q\ = %\Ed\- Therefore, 
for the current to be stable \Ed\ < A. This is the 
Landau criterion of superfluidity. For a tight binding 
lattice dispersion relation, we rewrite the condition as 
J sm(qa) smkpa < A. To complete a Bloch oscillation, 
sin(ga) should achieve its maximum value 1, i.e. 

sin Ufcl < A/ J. (6) 

For weak coupling, A/ J is given by the BCS theory, and 
in the attractive Hubbard model in the strong coupling 
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limit the gap at T = is given by A = tj- for half filling 
[li) . Using these estimates, we show in Fig.|5]the relation 
© for a gas of 6 Li atoms together with the transition 
temperature. To relate the criterion to the Cooper pair 
size, we rewrite Eq.© in terms of the BCS coherence 
length £0 = and insert Jsin(fcpa) = %vp/a which 
yields £0 < a/ 71 "- The observation of Bloch oscillations is 
thus restricted to superfiuids with BCS coherence length 
smaller than the lattice periodicity. This is the inter- 
mediate - strong coupling regime. The length argument 
can be also understood by thinking that the pairs have 
to be smaller than the lattice sites in order to see it as a 
periodic potential. 

For calculating the superfluid velocity a space depen- 
dent description of the superfluid has to be used. We 
combine the BCS ansatz with the Bloch ansatz for the 
lattice potential using the Bogoliubov - de Gennes (BdG) 
equations [T^j- As given by the Landau criterion above, 
the interesting regime is the intermediate - strong cou- 
pling one. Note that even in the strong-coupling limit, 
the algebra of the BCS theory can be applied to all 
coupling strengths 0, 0] together with an extra defi- 
nition for the chemical potential which in the weak cou- 
pling limit is given just by the Fermi energy of the non- 
interacting gas. The BdG equations are: 




FIG. 2: The transition temperature, Landau criterion at 
T=0 and the amplitude of the velocity Bloch oscillations for 
e 'Li atoms in hyperfine states with scattering length a s — 
-2.5-10 3 a for a half filled 3D CO2 laser lattice (a = 10 5 a )as 
a function of the lattice depth. The amplitudes of the oscilla- 
tions at T = 2/3T™ ax (horizontal line) are denoted by * for 
the normal state oscillations, o for the superfluid velocity at 
T=0 Eg. 1121 and in the boson limit Eg. 11131 pair size I — a/3 
by x and + for pair size / = a/4. The Landau criterion con- 
dition Eq.(|HJ requires A/J > 1 for the half filled band. Here 
Er is the recoil energy and E F is the Fermi energy for free 
fermions with the same density. 



H{v) - /x A(r) 
A(r)* -[H{t)-h] 



u(r) 
v(r) 



= E 



u(r) 
v(r) 



(7) 



When the external potential is periodic one can use the 
Bloch ansatz for u and v because by self-consistency the 
Hartree and pairing fields are also periodic. We obtain 



ttk(r) 



• r «k^k(r) ; tfc(r) = e 4k ' r « k k (r) (8) 



solutions for which -E^ q 



> 0, i.e. mm(y/E 2 A + |A q | 2 ) = 
|A q | > Ed- Remarkably, this condition is the same as 
obtained using the Landau criterion. 

In the BCS ansatz, a common momentum q can be 
added to all particles, leading to correlations of the 



type (c k+q c^ k+q ). The superfluid net momentum be- 
comes 2q. One can formally calculate this obvious result 



A 0) = }Z \9\ [1 - 2/(£ k )]w k (rK(r), (9) a i so by using the plane wave ansatz u k = |u k |e*( k+q )' r , 
k Vk = |u k |e i ( k+q )- r [13 (Eq.JTD} with </>=!) and intro- 



where k are the Bloch enveloping functions, such that 
[H(r) - A# k e 4k - r = a^ke 4k ' r . 

To describe Bloch oscillations we impose the adiabatic 
condition, that is, momenta evolve according to k — ► 
k + Fta/h = k + q, i.e. we consider BCS state with a 
drift (again only in x-direction) . The solutions of the 
BdG equations take the form 



A q (r) = 



— „*-r iq-p-q 



e e 



kVk+ q (r );« k (r) = e 



-iq-r~q 



k 



. 9 |[l-2/( J B+ q )] U q (rX(r)e 2iq - r , 



K = K 
Ed 



k+q 



-Jq-k)/2 ± V / (a +q +eq-k) 2 /4+|A q | 2 = 

where Ed is the energy differ- 



|A q | 



ducing an (unnormalized) order parameter wave func- 
tion A q (r) = e i2qr C, where C is given by Eq.iO to 
be a constant in r. Expectation values like momen- 
tum (p = —id/dr) can be calculated: (p) = (A q (r)| — 
z<9/<9r|A q (r))/(A q (r)|A q (r)) = 2q. The order parameter 
wave function is defined in the spirit of (but not with a 
one-to-one correspondence to) the Ginzburg-Landau the- 
ory with a space dependent wave function whose abso- 
0k- q (r) lute value equals the gap. In case of Fermionic atoms 
the Ginzburg-Landau approach has been used to describe 
harmonic confinement [l9| and vortices |20j . For the pe- 
riodic potential we introduce the order parameter wave 
function in the form |A q (r)) = ^ k |A q (r)), where using 



(10) 



ence and Ea the average energy. The ± holds for 
the particle and hole branch, respectively, and the 
particle branch eigenfunctions are |« q | 2 , |v k | 2 = (1 ± 
E A /\/E A + |A q | 2 )/2. The Hamiltonian transformed un- 
der the Bogoliubov transformation leading to (JJJ has to 
be positive definite. This means that one should use the 



|A q (r)) =^(k,q)|0 k+qe l ( k + q ) r )|<_ qe -( k - q >) (11) 



and F(k,q) = \g\[l - 2/(£ k hq )]u q 6 q * . We calculate 
the superfluid velocity using (v s ) = A/"(A q (r)|r|A q (r)), 
where N = (A q (r)|A q (r))- 1 . Using (r) t e -, (k - q) . r = 
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IcHc^k-q an( i *he tight-binding energy dispersion rela- 



tion the superfluid velocity becomes 
(vxs) = A/~y^ |F(k, q)| 2 ^cos k x a sin qa 



^sin (qa)Af^2 



[1 - 2f{E*)\ A q 



I Ail 



(12) 



cos fc x a. 



The superfluid velocity for selected parameters is shown 
in Fig. [21 We have also calculated the thermal quasipar- 
ticle contribution but is turns out to be negligible. 

In the composite boson limit, one could describe the 
center-of-mass movement of the composite particle by 
defining J* = J(m — » 2m). In order to give a simple 
estimate for the effect of the Fermi statistics, we inter- 
pret |F(k,q)| 2 ~ \F(k)\ 2 in Eq.^ as reflecting the in- 
ternal wavefunction of the pair in the composite boson 
limit, c.f. 0, 0]. The average velocity for the bosons 
becomes (v x b) oc sinqa^ k |i^(k)| 2 cosk x a. If the 
pairs were extremely strongly bound, the internal wave 
function in real space is a delta-function, corresponding 
to a constant in k-space. This means (vb) — since 
the cosine integration in Eq. l|12l) would extend to the 
whole k-space with equal weight, i.e. there are no empty 
states in the Brillouin zone as required for Bloch oscil- 
lations. For on-site pairs, we use l-F^)! 2 oc exp(— r 2 /l 2 ) 
leading to \F(k)\ 2 oc AT exp(-l 2 k 2 /4), therefore the sup- 
pression factor for the Bloch oscillations becomes S ~ 
Af J dk exp(— l 2 k 2 /A) cos ka, where I is the pair size. As 
a rough estimate for the average velocity we thus obtain 



sin (F x ta/K). 



(13) 



This is shown in Fig. [21 for pair sizes I = a/3 and I = a/4. 
It gives an order-of-magnitude estimate, approaching the 
results given by the BCS algebra. 

Another way of treating the composite boson limit is 
to derive a Gross-Pitaevskii type of equation for the com- 
posite bosons with M — 2m and with a repulsive non- 
linear interaction term ubUb = ns^irf^aB/M, as = 2a s 
where a s is the renormalized s-wave scattering length 
[2l| . If the non- linear term is small compared to the 
Bloch energy Eb = h 2 /(Ma 2 ), the nonlinearity leads 
only to a change in the band width J |3j. Therefore, com- 
posite bosons oscillate but with a modified amplitude. 
Large non-linearity would not allow Bloch oscillations, 
corresponding to a large suppression factor in the above 
discussion. Note that the Landau criterion for a super- 
fluid Bose gas gives the critical velocity v soun d - 



M 



which is orders of magnitude bigger than Eq. (| 1 31) for half 
filling and parameters in Fig. |2 Problems arise only in 
the extremely empty lattice limit. 

In summary, we have defined a set of tools for quali- 
tative and quantitative description of Bloch oscillations 
for the BCS-BEC crossover regime. The amplitude of 



the oscillations decreases when the crossover is scanned, 
in general due to the shrinking of the bandwidth. How- 
ever, the change from the normal to the superfluid state 
description leads to a drastic change in the amplitude. 
This is due to smoothening of the Fermi-edge by pair- 
ing. Bloch oscillations could be used for exploring pair- 
ing correlations since any localization in space (pair size) 
leads to broadening in momentum which suppresses the 
amplitude in the same way as band filling in the non- 
interacting gas. Achievement of superfluidity is still a 
great challenge, but even at T » T c , the effect of 
collisions on Bloch oscillations can be studied produc- 
ing information useful for applications of Bloch oscilla- 
tions such as production of Terahertz radiation [lfll l22j | . 
Observation of oscillating fermionic atoms in optical lat- 
tices would contribute to the quest for a steadily driven 
fermionic Bloch oscillator. 
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